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МОТОРНОГО ТОПЛИВА В ДИЗЕЛЬНОМ ДВИГАТЕЛЕ" 
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(Донской государственный технический университет) 


В целях выбора адекватного описания процессов впрыска, распыливания и испарения моторных топлив в 
дизелях предложена иерархия соответствующих МОДЕЛЕЙ. Исследована область применения детальной мо- 
Дели, основанной на эволюционных уравнениях с частными производными, а также осредненных вариантов, 
выражаемых посредством двух и одного обыкновенных дифференциальных уравнений. Обосновано исполь- 
зование упрощенных МОДЕЛЕЙ в компьютерной имитации рабочем О процесса ДИЗЕЛЕЙ с разрешенным описа- 
нием образования ТОКСИЧНЫХ газов И сажи, 
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Введение. Вклад автотранспорта в загрязнение воздушной среды мегаполисов превышает 90 %. 
Кроме того, рост российского автопарка примерно на 60 % обусловлен его старением. Отсюда 
и потребность оценивать энергоэкологические характеристики двигателей внутреннего сгорания 
(ДВС) в зависимости от технического состояния, в частности состояния топливной системы, а 
также состава моторного топлива [1 — 3]. Реалистичное компьютерное моделирование выступает 
действенным инструментом, позволяющим совершенствовать существующие и создаваемые 
транспортные ДВС. 

Постановка задачи и описание моделей. Способы формирования заряда в ДВС в значитель- 
ной степени определяются агрегатным состоянием и другими физико-химическими свойствами 
топлива. Если в качестве энергоносителя используется каменный уголь, требуется произвести 
предварительную специальную обработку и подготовить топливо к сжиганию: измельчить его до 
пылеобразного состояния и смешать с жидким горючим или транспортной средой (водой) в 
эмульсию (суспензию). Двухфазную систему — однородную взвесь призмовидных твердых частиц 
углерода в жидкости [4] —- впрыскивают в камеру сгорания (КС) через форсунки. Таким образом, 
технология сжигания в ДВС твердых угольных топлив представляет собой модификацию рабочего 
процесса с внутренним смесеобразованием. 

При внешнем смесеобразовании в КС попадает сравнительно однородная по химическому 
и фазовому составу топливно-воздушная смесь, что позволяет квалифицировать динамику такого 
процесса как вырожденный случай внутреннего смесеобразования -— с бесконечно быстрыми фа- 
зовыми переходами и гомогенизацией в начале такта сжатия. Поэтому, выбрав в качестве базо- 
вой модель внутрикамерного смесеобразования с твердыми частицами в составе топлива, мы 
предусматриваем все практически реализуемые ситуации. 

Подача топлива в цилиндр осуществляется из трубопровода высокого давления Р,, а ее 
динамика определяется изменением во времени проходного сечения 9. Последнее служит пара- 
метром управления процессами в ДВС с внутренним смесеобразованием. Массовая скорость 
впрыска составляет [5]: 

КО = 9(0) [2(Р.-Р)/р]" = (2Р/р)"? (В, (1) 
где Р- давление рабочего тела в КС; р — плотность. 





* Работа выполнена при финансовой поддержке фонда АЁСОА. 
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В элементах форсунки поток турбулизуется и, истекая в заполненную газом КС, дробится 
на капли вследствие тейлоровской неустойчивости [6, 7]. Распаду жидкой струи на капли препят- 
ствует энергия поверхностного натяжения Гельмгольца. Отношение мощности формовозмущаю- 
щего и стабилизирующего факторов известно как число Вебера: 

\!е = ри’ Л/оцьь, (2) 
где и - скорость капель; Л - диаметр капли; оз» — коэффициент поверхностного натяжения. 

Установившемуся дисперсному составу капель отвечает диапазон чисел \е = 9-18, т. е. 
на диспергирование струи тратится 5-10 % взаимной кинетической энергии фаз [5, 6]. Возмуще- 
ние жидкости в системе топливоподачи и турбулизация газа в КС направлены на улучшение каче- 
ства распыла, т.е. на возможно более полное и мелкое дробление струи. Согласно [8], нагрев то- 
плива сопровождается уменьшением коэффициента поверхностного натяжения, что также улуч- 
шает распыл: 

Опов — а (1 И Вр: (3) 
где Т„»- критическая температура; 0’в- экстраполированное значение опь(Т — 0). 

Критическое значение числа Вебера определяется в эксперименте с точностью до множи- 
теля - 2, что обусловлено полидисперсностью распыла. Шаблоном распределения капель по разме- 
рам выберем аналитически интегрируемый вариант обобщенной зависимости Розина-Раммлера [9]: 

Ерь(А) = а Л ехр(-УЛ°), (4) 
где а, В, уи д - экспериментально определяемые константы. 

Подобрать три параметра модели (4) по экспериментальным данным, объем и точность 
которых ограничены, проблематично. Требование аналитической интегрируемости (4) позволяет 
связывать В и б соотношением В = 6 - 1. Неопределенность критического числа Вебера составляет 
=2!? [5, 6]. Распределение (4) имеет также подходящую ширину по уровню 0,1'тах { Ер-ь(А)} при 
б = 3. Наконец, распределение (4) с 6 = 3 почти симметрично относительно максимально вероят- 
ного значения и близко к нормальному. 

Предложенная модель дробления струи позволяет связать показатель диспергирования 
Еъ.р(А) с единственным управляющим параметром: 

у = У(<\е>) = 2/3 [ри’Конь<Ме>)[ = 3:10“ [ри’/онеь], (5) 
который следует трактовать как обратный характерный размер структурно устойчивых частиц 
жидкой фазы. 

Из уравнений (1) и (2) следует, что начало и конец впрыска характеризуются наихудшим 
измельчением струи. Практически крупные капли раннего формирования успевают испариться и 
сгореть, а финальная порция плохо распыленной жидкости оказывается в невыгодных условиях. 
Отмеченное обстоятельство определяет неполноту сгорания и оправдывает полезность следую- 
щей оценки. Для модели монодисперсного распыла мгновенный диаметр капель 


А -и?-9(0”. (6) 
Когда пропускное сечение открывается мгновенно в момент 1.» 
9(2) тя $121(Е ы Ь (7) 


проблем не возникает: все капли идентичны. В противном случае часть струи измельчается не- 
достаточно. Долю неудовлетворительного распыла можно определить, приблизив начальный и 
конечный участки функций Р(@ и 9(1) степенной функцией вида Г: при несимметричной динамике 
впрыска параметр «+, конечно, отличен для этих участков. Диаметр структурно устойчивых жидких 
частиц, образующихся впоследствии из порции потока, впрыснутого в промежутке времени 
#1 (#+ 40), Л(® пропорционален #2“, а их общая масса — 3() 4. Обозначив максимально допусти- 
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мый диаметр капли Л„‚, путем несложных преобразований получим зависимость массовой доли 


недораспыленного топлива: 

ТЯ — (АзиЛА а) НВ 5 (8) 
где Л — диаметр капель в момент максимально интенсивного впрыска &„„„. При выводе (8) пред- 
полагалось, что впрыск с различной динамикой (значениями о) характеризуется одинаковыми 
продолжительностью („„„„ и суммарной дозой (массой) жидкого топлива, причем Л,, > Ли». Техни- 
ческая подоплека этих ограничений очевидна: а) для эффективного смесеобразования требуется 
осуществлять его в течение ограниченного периода; 6) если чрезмерно крупными оказываются 
даже капли, отвечающие наилучшим условиям дробления, повышение качества распыла требует 
наращивания энергетики впрыска. 

Интеграция моделей впрыска и неоднородного диспергирования сопровождается появле- 
нием «перекрестных членов», отвечающих мультипликации «слабых мест» обоих процессов. По- 
этому динамику впрыска следует оптимизировать очень аккуратно. 

Распыливание топлива в КС сопровождается постепенным сопряжением траекторий от- 
дельных капель с траекториями мезообъемов газов. После измельчения дроблением ведущая 
роль отводится процессам испарения жидких частиц и газофазной гомогенизации и/или выгора- 
ния отдельных капель. 

Бесполезность наиболее общей модели испарения жидкого топлива для наших целей ста- 
новится очевидной, если учесть следующие обстоятельства [4, 10, 11]: 

— испарение может происходить в до- и закритическом режимах; 

— многокомпонентные топлива могут испаряться как с «разгонкой» фракций, так и без 
нее; 

— испарение может сопровождаться крекингом топлива или происходить без реакций; 

— фронт горения может примыкать к поверхности испарения или отсутствовать в диффу- 
зионной окрестности; 

— возможны ламинарный и турбулентный режимы обтекания жидких частиц; 

— в состав жидкого топлива могут входить мелкодисперсные частицы твердых горючих 
веществ, например угольная пыль; 

— эмульсионные топлива часто содержат воду, критическое давление которой на порядок 
выше, чем Р‚„, используемых углеводородных топлив и выше достижимого давления в КС. 

Корректно учесть перечисленные обстоятельства не представляется возможным, поэтому 
динамику фазовых превращений будем описывать с опорой на особенности испарения жидкого 
топлива в типичном ДВС, работающем по циклу Дизеля [10, 11]. 

Б. Срезневским экспериментально установлено, что диаметр испаряющейся сферической 
капли Л уменьшается по закону: 

а =- р. А", (9) 
где Ш, - коэффициент диффузии паров топлива в вакуум. 

Коэффициент диффузии зависит от молекулярного веса топлива ц., температуры и давле- 
ния газа в КС: 

О” Е (10) 
если давление измерять в атмосферах, а Р» - [см?/<]. 

Динамика газификации топлива, поступающего в КС с массовой скоростью С(1), задается 
уравнением: 

90—12 А 9ЙОЛ = ЕСА, 8, (11) 
которое связывает темп фазового перехода с дисперсностью жидкости (ее определяет функция 
распределения капель по размерам /= ХА, 0), мкм ') и законом топливоподачи РКА, 2), мкм"*.с". 

Для наших целей представляет интерес связь между массой испаренного на некоторый 
момент топлива т(р) и динамикой его впрыска С(0. Последовательное рассмотрение вопроса 
включает два этапа: первый — решение уравнения (11) для конкретного источника РКА, 0}; второй — 
вычисление мгновенной скорости накопления жидкого топлива в КС и темпа газификации капель 
на основе КА, 41). 
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Семейство характеристик (11) задается уравнением: 
ал =- р. /А или Л^=Л-Ь:Ь (12) 
гДе Ло отвечает начальному размеру капли. Действие источника ЕКА, 2) сводится к «загрузке» ка- 
плями соответствующих характеристических линий. Сингулярность - 1/Л в коэффициенте уравне- 
ния (11) при ОДОЛ побуждает прибегнуть к замене независимой переменной Л —> т = р? и опе- 
рировать уравнением: 


9101 — 3/2 р, т дЙдт = М(т, 9. (13) 
Соответственно изменятся функция распределения капель по размерам 
Ер (Л) = ЗУЛ* ехр(-УЛ°) — — Мьыт)= р ехрС-/рт) (14) 
и уравнение характеристик (ниже константа Л, увеличена в 6 раз): 
ЕЮ 55 моет” = Вереи (15) 


Нормированная на единицу величина Мь (и?) во втором из уравнений (14) легко модифицируется 
для учета временного фактора. 
Учитывая, что согласно (4)—(6) показатель измельчения струи у — С(0)°, получаем явный 
вид источникового члена в уравнении (13): 
Мп, 0 = (/РУ’С(" ехр(-Е/р С(0° т), (16) 
где &- коэффициент, характеризующий распыливающее устройство. 
Нетрудно убедиться, проинтегрировав (16) последовательно по обеим переменным, что 


| т М(т, ) т=с@ и | С(® &= (Ц;. (17) 
0 0 
В интегралах (17) нормирование осуществлено по цикловой подаче — массе впрыскивае- 
мого за рабочий цикл топлива С». 
При выбранных модельных распределениях средняя масса капель <т({)>, поступающих в 
КС в момент времени 1: 
<т(0> = р/и=р/[5 (0, (18) 
что будет использовано при упрощении модели впрыска-распыливания. 
Уравнения (11) и (13) позволяют физически корректно учесть динамику испарения капель 
неодинакового размера, поступающих в КС в различные моменты времени. Общее решение (13) 
получается интегрированием функции источника Мь р(т, #) вдоль характеристик (15) и имеет вид: 


1 
Кт, = | Мьь ([т?3 + Ш, р?3&- ©.) а. (19) 
0 
На основании массово-геометрических характеристик жидкой фазы в КС (19) можно опре- 


1 
делять массу испаренного топлива как разность масс впрыснутой [6(%т и еще не испаренной 
0 


[2] 


| Кт, 8) ат жидкости. Рассчитываемая таким способом массовая скорость газификации, состав- 
0 


ляет: 


Го = | т[дКт, В/О ат. (20) 
0 
Упростить модель позволяют следующие соображения. Правую часть (20) можно записать 
в виде: 


| т[ЭКт, 2/91] ат = д/0 | т т, 9) ат = 9/9 {<т()> | Кт, ат. (21) 
0 0 0 
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Если дополнительно предположить, что характерное время испарения отдельных капель 
меньше, чем временной масштаб изменения их числа, можно представить правую часть (21) как 
о о 
9/0 {<т(р> | Кт, в ат} = д<т(И>/0Е | От, в ат. (22) 
0 0 

Оценить первый множитель в (22) можно при помощи уравнения характеристик (15): 
д<т(®>/ди = 3/2 ` р, р’? <т(®> '3. Второй, представляющий собой число капель в КС, очевидно, 
равен М(0/<т(#)>, где М (1) — масса жидкой фазы. 

С учетом выполненных замен и подстановок (20) трансформируется в уравнение: 

Г = 3/2.) (р/<т(®>)* МО, (23) 

а динамика суммарной массы жидкой фазы в КС определяется уравнением (23) и условием мате- 
риальной консервативности: 

аАМ/ = С(а) — 3/2. (р/<т(®>)? М. (24) 

Подобная модель с постоянной величиной <т(#)> = <т> использована в работе [12]: по- 
лидисперсность распыла учитывалась качественно — заменой «сингулярного» уравнения испаре- 
ния, отвечающего конечному времени существования капли, «размытым», согласно которому ка- 
пли существуют вечно. 

Далее следует задать правдоподобный закон эволюции параметра (<т(#)>/р)?°, отвечаю- 
щего площади поверхности среднестатистической капли в КС: необходимо учесть влияние на нее 
как собственно испарения, так и влияние источника. Нетрудно убедиться, что уравнение 

ф<т(о>” “1 = — Бр ^® + [р *°/8 *°6(0*) — <т(0 > Е СОРМ то), (25) 
в котором размер вновь поступающих капель определяется законом впрыска С(1) и соотношением 
(18), удовлетворяет предъявленным требованиям. 

Система обыкновенных дифференциальных уравнений (ОДУ) (24)—(25) представляет со- 
бой упрощенную модель испарения топлива в дизелях. Для ее получения использовался метод 
осреднения уравнения с частными производными (УЧП), состоящий в понижении пространствен- 
ной размерности путем введения «эквивалентной» нелинейности. 

Ценой отказа от использования в модели испарения УЧП (13) служат жесткость и нели- 
нейность результирующих ОДУ (24)-—(25), заметно усложняющие их численное интегрирование. 
Кроме того, в момент начала впрыска и/или при чрезвычайно интенсивном испарении М = 0 ди- 
намика средней площади капли не подчиняется (25) и необходимо пользоваться уравнением (18). 

Согласно (24), в динамике жидкой фазы выделяются переходный и квазистационарный 
этапы: в начале подачи топлива М (9 = | С(®) Е; в установившемся режиме М(®, подстраиваясь 
под темп впрыска, принимает значение около 2/3, р 2“ ?С(1*]; по окончании поступления 
капель в КС величина М стремится к нулю с характерным временем РБ; '(<т(0>/ру". С другой 
стороны, такое поведение обеспечивается свойствами уравнения (13), характеристики которого 
разбивают мировую плоскость {т х 1 на области влияния источника и начальных данных. Пере- 
ходные режимы в динамике жидкой массы, определяемой ОДУ (24)-(25), обусловлены простран- 
ственно-временной ограниченностью источника капель. 

Функцию впрыска С(1) удобно задать линейной комбинацией У, С{1) звеньев вида: 

СК = А, (о, + 1) (1- В) АЕ" ехр[-@ - Ав)" вп — 6), (26) 
где А, — удельный вес; а, №; и Дё — форм-факторы 1-го звена. 

Сумма нормируется очевидным образом: Х 4; = С;. Практически для аппроксимации ре- 
альных законов топливоподачи достаточно двух слагаемых типа (26). Удобно положить о = 1/3, 
02 = 3, = идя АН, аА,/А› выбрать в пределах 3-5. Если имеются данные о реальной дина- 
мике впрыска, целесообразно параметры модельного распределения (26) определить как пара- 
метры аппроксимации. 
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Результаты компьютерного эксперимента. Для сопоставления различных способов описания 
кинетики газификации жидкого топлива в дизелях приведем результаты соответствующих тесто- 
вых расчетов. Объектом моделирования выбран цилиндр тепловозного дизеля: его рабочий объ- 
ем составляет (15 - 20) . 103 смз, частота вращения коленвала - 10 - 15 с ', цикловая подача топ- 
лива С’. — 1 г, ДЁ-- 4 мс, характерный начальный диаметр капель топлива -— 50 мкм, коэффициент 
диффузии 2» — 10? см?/с. 

Отвечающая выбранным параметрам динамика впрыска и распыливания капель представ- 
лена на рис. 1. 
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Рис. 1. Количественные характеристики впрыска топлива в тепловозном дизеле: 
а - массовая мгновенная подача жидкости; 6 -— массовая интегральная подача жидкости; 
в - средняя масса капель; г - секундное число капель; С' = 0,77 ги С. = 0,23 г - массовые доли 
топлива, относимые соответственно к слагаемым формулы (26) 


Данные рис. 1 свидетельствуют, что качество смесеобразования определяется продолжи- 
тельностью и переходными режимами топливоподачи. 

Теперь сопоставим описанные модели распыливания-испарения различной детальности с 
тем, чтобы выбрать адекватную для интеграции в новую версию автоматизированной системы 
МОТОР [13, 14]. 

Наиболее подробно процесс газификации топлива описывается уравнением (13). Его ре- 
шением (19) для выбранных нами условий, служат массово-временное распределение частиц в 
источнике, камере сгорания, а также момент т т, #), изображенные на рис. 2. Условия расчетов 
соответствуют тепловозному дизелю на режиме номинальной мощности. 
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Рис. 2. Статистические характеристики распыленного топлива в дизеле согласно модели (13), (16): 
а- масса жидкого топлива в камере сгорания; 6 -— распределение капель по размеру; 
в- «массовая представительность» капель размера т 


Результаты компьютерного эксперимента, показанные на рис. 3 и 4, призваны прояснить 
зависимость динамики испарения топлива от качества распыла, задаваемого величиной &, и ко- 
эффициента диффузии топлива Д.. На рис. 3 видна геометрия зон влияния этих параметров на 
конечный результат. В частности, отражены: 1) усиление зависимости Дт, д) от О, на завершаю- 
щей стадии процесса (рис. 3, а); 2) устойчивость формы распределения капель по размерам при 
разумном варьировании показателя диспергирования &; 3) область сильного влияния одновре- 


менно обоих параметров (характеризуется высокими значениями и сгущением изолиний) совпа- 
дает с завершающей фазой впрыска. 
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Рис. 3. Влияние коэффициентов р, (а) и & (6) на распределение капель 
по размерам согласно модели (13) в условиях рис. 1и 2 
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Рис. 4. Влияние коэффициентов р, (а) и & (6) 
на динамику «массовой представительности» капель 
согласно модели (13) в условиях рис. 1и 2 


Для ответа на вопрос, насколько правомочным является упрощенное описание впрыска 
при помощи ОДУ (23)-(25) и их предельных форм, необходимо выполнить специальные сопоста- 
вительные расчеты. В качестве сопоставляемых моделей выберем следующие три: 1) основанную 
на УЧП (13); 2) использующую аппроксимирующие ОДУ (23)-(25); 3) совпадающую с предыду- 
щей, но предполагающую постоянство средней массы капель, и в результате определяемую един- 
ственным уравнением (24). 

Естественными объектами такого сопоставления (его результаты представлены на рис. 5 и 6) 
являются усредненные по мгновенному массовому распределению К, В) величины <т(й>?? и 
М (8. На рис. 5 показаны средняя поверхность испаряющихся капель (рис. 5,а) и масса жидкого 
топлива (рис. 5,6) в КС тепловозного дизеля, рассчитанные в рамках различных моделей испаре- 
ния: сплошные линии без маркеров — модель, основанная на УЧП (13): заштрихованный на фраг- 


менте рис. 5,а коридор соответствует интервалу <т?/р> — <т/р>??. 











(тбдЬ/ р)" 105 см? М, 
в - 0,3- 
4 - 0,2- 
я 0,1- 
0 5 101, мкс 0 5 10 #, мкс 
а) 6) 


Рис. 5. Характеристики распыливания-испарения согласно (13), (24)-(25), (24): 
а- средняя поверхность капель <т(0>?3; б- количество жидкого топлива в камере сгорания; 
х — описание при помощи системы ОДУ (24)-(25); П — модель, использующая одно ОДУ (24) 
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Представленные на рис. 5 данные подтверждают квазистационарность среднего размера 
капель в КС, чем оправдывают использование предложенных здесь упрощенных моделей испаре- 
ния. Погрешность упрощенных моделей газификации топлива (24) и (24)-(25) в редких случаях 
превышает 20 - 30 %, что позволяет использовать в интегрированной модели рабочего процесса 
ДВС математическое описание испарения, основанное на предположении о постоянстве среднего 
размера капли и выражаемое единственным ОДУ (24). 

Существенно, что сопоставляемые модели, несмотря на различную степень детальности, 

сходным образом воспроизводят зависимость доли испаренного (за отведенное время) топлива от 
качества распыливания & и легкости испарения Л. Расчеты показывают, что неполнота испаре- 
ния, а следовательно, и сгорания сильно зависят как от качества распыливания (ОшМ.„Ошё —> —1), 
так и от транспортных свойств рабочего тела (ОшМ»..„О\ Г —> 2). 
Выводы. Совпадение динамики испарения в предложенной иерархии моделей различной под- 
робности заметно облегчает моделирование внутрицилиндровых процессов, в том числе образо- 
вания токсичных веществ. Результаты работы служат также строгим научным обоснованием уп- 
рощенного описания испарения посредством модели (24), применявшегося авторами в предшест- 
вующих исследованиях. На основании предложенных моделей выявлены критичные условия и 
режимы в системе процессов впрыска, распыливания и испарения топлива в дизелях. 
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